Staggered pairing phenomenology for UPd 2 Al 3 and UNi 2 Al 3 



V. Martisovits 

Department of Physics, Ohio State University, Columbus, Ohio 43210 

D. L. Cox 

Department of Physics, Ohio State University, Columbus, Ohio 43210 r-, 
and Department of Physics, University of California, Davis, California 956 laJ 

We apply the staggered-pairirig Ginzburg-Landau phenomenology to describe superconductivity in 
UPd2Al3 and UM2AI3. The phenomenology was applied successfully to UPt3 so it explains why 
these materials have qualitatively different superconducting phase diagrams although they have the 
same point-group symmetry. UPd2Al3 and UM2AI3 have a two-component superconducting order 
parameter transforming as an ff-point irreducible representation of the space group. Staggered 
superconductivity can induce charge-density waves characterized by new Bragg peaks suggesting 
experimental tests of the phenomenology. 



Heavy-fermion superconductorsBi exhibit many signs 
of unconventional superconductivity thus playing an im- 
portant role in advancing our knowledge of this impor- 
tant subject. One of the clearest signs of this unconven- 
tional pairing is a complex phase diagram with multiple 
superconducting phases. This is clearly the case for IJPt3 
which has three distinct superconducting phases uu On 
the other hand, UPd 2 Ai3 and UJ^ALj have simplft-phase 
diagrams with only one superconducting phaseJa The 
crystal structures of UPt3, UPd2Al3, and UM2AI3 are 
characterized by the same point group (Dq^) and the 
same Bravais lattice (simple hexagonal) with different 
numbers of uranium atoms in the basis (two for UPt3, 
one for UPd 2 Al 3 and UNi 2 Al 3 ). 

So far no comprehensive microscopic theory exists_fpx 
these materials. Many phenomenological approachesLTtHl 
have tried to explain exotic superconductivity in UPt3. 
Most of them use a multidimensional superconducting 
order parameter transforming as a representation of the 
point group Dq^ and a Ginzburg-Landau free energy in- 
variant under this group. If one tries to apply this ap- 
proach to describe superconductivity in UPd2Al3 and 
UM2AI3, it is necessary to choose an order parameter 
transforming as a different representation of the group 
Dgh- No single phenomenology is able to predict satis- 
factorily the relevant representation for each material. 

In this paper, we 

apply the staggered-superconductivity Ginzburg-Landau 
theorylljiij to describe superconductivity in UPd2Ai3 and 
UNi 2 Al3. This theory has been successful in explaining 
the H-P-T superconducting phase diagram of UPtsETEl 
including the tetracritical point for all directions of the 
magnetic field. The distinctive feature of this theory is a 
superconducting order parameter transforming as an irre- 
ducible representation of the space group characterizing 
the crystal structure which resides away from the center 
of the Brillouin zone. Namely, we have finite center-of- 
mass or staggered pairing. For UPd2Al3 and UM2AI3, 
the star in the reciprocal space carrying the representa- 
tion is determined by the negative pair hopping model. 



We introduce a two-component -ff-point superconducting 
order parameter and write a free energy including gradi- 
ent terms and coupling to a symmetry-breaking strain. 
Minimizing the free energy in the absence of the magnetic 
field, we find stability domains in the parameter space, 
and identify sequences of possible phases for each domain. 
Application of an external or spontaneous symmetry- 
breaking strain does not split the superconducting tran- 
sition. The calculated upper critical field as a function 
of the temperature shows no kink. We introduce charge- 
density- wave and magnetic order parameters and investi- 
gate their possible coupling to the superconducting order 
parameter. Charge-density-wave order can be induced by 
staggered superconductivity suggesting a possible exper- 
imental test of the phenomenmology. On the other hand, 
magnetic order cannot be induced. 

The motivation for the staggered-superconductivity 
theory is based upon the theory of odd-in-frequency pair- 
ing characterized by a gap.riunction which is an odd 
function of the frequency.t2lli-3 This motivation is two- 
fold: (i) there are strong reasons to believe that an odd- 
in-frequency superconducting state cannot .be stable for 
pairs with zero center-of-mass momentum and (ii) 
the two-channel Kondo lattice, which is one candidate 
theory for heavy fermion compounds, has recently been 
shown to produce odd-frequency superconductivity in in- 
finite spatial dimensionalJ which is likely to be staggered 
for finite dimensionality. (For a recent review of tato- 
channel Kondo lattice models and related physics, seeEj.) 

In consequence, we develop a phenomenology for a 
non-uniform superconducting state described by a stag- 
gered superconducting order parameter ancL-a negative 
pair hopping energy between uranium sites. 113 We then 
write a hopping model for Cooper pairs assuming in- 
plane and out-of-plane nearest neighbors. The energy 
band for this model has a two-fold degenerate minimum 
for the crystal momenta located at the opposite corners 
H\, H2 of the Brillouin zone shown in Fig. |. (For UPt 3 , 
the different number of uranium_atoms in the basis moves 
the minimum to the M points.til) Cooper pairs will con- 
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dense in states with these momenta so in this picture 
there are pairs with finite center-of-mass momenta. The 
Bloch functions ^i(r),^ 2 (r) of these states are uniquely 
determined by specifying the Wannier function. We as- 
sume that it transforms as a real one-dimensional rep- 
resentation of the uranium-site point group D§h- For 
example, it will transform as A2 if we use a two-.channcl 
quadrupolar Kondo model to describe UPd2Al3.E£l How- 
ever, all the following results do not depend on this spe- 
cific choice. 

We introduce a two-component superconducting order 
parameter 771,772 using the gap function A(r)t2l where r 
is the pair center-of-mass coordinate: 



A(r) = (771 - ir) 2 )ipi(r) + (771 + ir) 2 )ip 2 (r). 



(1) 



The choice of the coefficients in front of the Bloch func- 
tions gives the usual transformation properties under the 
time reversal transformation 771 — > 77^ , 772 — > . With re- 
spect to the space group Dh, the order parameter trans- 
forms as a two-dimensional i7-point irreducible represen- 
tation. 

The free-energy density F is required to be invari- 
ant under the elements of the space group D\ h , the 
U(l) gauge transformations, and time reversal. We 
include coupling to the hexagonal-symmetry-breaking 
strains e\ ~ x 2 — y 2 , e 2 ~ 2xy which transform as the 
center-of-Brillouin-zone representation Ei g . The bulk 
part Fq and the gradient part F g are given by 



magnetic order at elevated temperatures which, at least 
for UP^ALj is not of the small moment variety found in 
UPt 3 ).H 

Minimizing the free energy, we can investigate possible 
superconducting phases for T < T c in the absence of the 
magnetic field. There are two regions of stability in the 
parameter space Pi , P2 in which the free energy has a 
stable lower bound. For the region I defined as j3\ > 0, 
P2 > there is only one phase (771,772) = (i, 1) breaking 
the symmetry under the time reversal (as favored by the 
last term in the bulk part F ). The region II defined 
as Pi + (i 2 > 0, f3 2 < is characterized by the phase 
(771,772) = (0, 1) having-pnly one component of the order 
parameter condensed. E2l 

Let us now discuss the upper critical field H C 2- Min- 
imizing the free-energy functional with respect to the 
components of the superconducting order parameter and 
keeping only terms linear in them, we get two identical 
decoupled linearized Ginzburg-Landau equations. This 
means that the field H C 2 does not have any kink as a 
function of the temperature T as observed experimen- 
tally (for UPt3 there are three different decoupled equa- 
tions, and the field H c2 is given by the largest of the 
three values generated by each equationrsp we get a kink 
corresponding to the tetracritical pointEJ). The depen- 
dence of H c2 on the direction of the magnetic field has 
the same form as the result for a conventional free energy 
with effective-mass anisotropy: 



F = F + F„ 



(2a) 
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F = a(r)(|77i| 2 + |772| 2 )+/3i(|77i| 2 + |77 2 | 2 ) ; 
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(2b) 



F 9 = K Y1 (\PxVn\ 2 + \PyVn\ 2 ) 

71 

+ K e £i ^ (\P^Vn\ 2 - \PyVn\ 2 ) 

71 

+ K e £2 ^ [(PxVn)(Pyrin)* + (PxVn)* (PyVn)} 
n 

+ KzJ2\P* T ln\ 2 - (2C) 



Here a(T) = a (T - T c ) with T c = T c0 - (e 2 + e 2 ) /a ; 
Pu P2, k 7 K ei and k z are phenomenological Ginzburg- 
Landau parameters; p a = —id a — 2eA a /Tic (e < 0). |?7i| 2 
and 1 772 1 2 enter the quadratic term in the bulk part Fq 
with the same coefficient a(T) which is valid up to any 
order in the strains £1,62 (the above expression demon- 
strating this explicitly up to the second order) . This im- 
plies that the presence of any finite symmetry-breaking 
strains do not split the superconducting transition which 
is in agreement with experiment in that only one super- 
conducting transition is found (despite the presence of 



R = (k 2 - K 2 |e| 2 ) cos 2 6 + k z [k- K £ Re (ee i2{p )] sin 2 9, 

(3b) 

with e — e\ — ie2- 

Finally, let us investigate the possibility of charge- 
density waves and magnetic order induced by staggered 
superconductivity. Starting from the staggered gap func- 
tion (1), we can introduce order parameters for them in 
a very general way.cll The gap function A(r) does not 
have the periodicity of the original simple hexagonal Bra- 
vais lattice; it has the periodicity of a new larger simple 
hexagonal Bravais lattice of Fig. ||with a' — \/3a, c' = 2c. 
Thus the most general charge and magnetic structures 
are described by this new lattice and a basis consisting of 
six uranium atoms. Assigning different charges and mag- 
netic moments for each element of the basis, wc define 
six-component order parameters. The six components 
can be replaced by three interplanar-bonding combina- 
tions and three interplanar-antibonding ones. To sim- 
plify the following discussion, we can ignore the bonding 
combinations because they do not couple to the super- 
conducting order parameter. Then the basis atoms in a 
plane and the basis atoms in a neighboring plane have the 
same charges pi,p2,P3 and the same magnetic moments 

Ml,M2,M3- 
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The components pi,p2,P3 of the charge order param- 
eter transform under the space group D^ h as a three- 
dimensional reducible representation composed of the 
totally symmetric irreducible representation and a two- 
dimensional if -point irreducible representation (Fig. ^). 
The totally symmetric combination p\ + p2 + P3 repre- 
sents the total charge in the large primitive cell which 
can be put equal to zero in the context of charge-density 
waves. Then the relevant part of the free-energy density 
is 

F c dw = a p (pj +pl + pt~ P1P2 - P\P?, - P2P3) 

+ g (\m\ 2 ~ \m\ 2 ) i^Pi - P2 - p 3 ) 

+ V3g (771772 + V1V2) (P2 - Ps.) (4) 

where a p and g are phenomenological Ginzburg-Landau 
parameters. Minimizing this free energy, we find possi- 
ble charge-density-wave structures for all the supercon- 
ducting phases identified above. There is no induced 
charge-density wave for the phase (771,772) = (*, 1) at the 
region I. For the phase (771,772) = (0,1) at the region 
II, staggered superconductivity induces a charge-density 
wave described by the order parameter {p\,p2,Pz) = 
(2, —1, —1). The structure is characterized by new Bragg 
peaks (|ni, ^774 + 77,2, 77,3) where 774, 77,2, and 77.3 are any 
integers. This important experimental prediction can be 
verified using neutron or x-ray scattering. 

The components fix , \i2 , A*3 °f the magnetic order pa- 
rameter transform under the space group as the prod- 
uct of the axial-vector representation and the represen- 
tation for the charge order parameter discussed above. 
The product is a nine-dimensional reducible representa- 
tion composed of two center-of-Brillouin-zone irreducible 
representations and two iGpoint irreducible representa- 
tions. In the free energy there are no invariants coupling 
the magnetic order parameter linearly to the square of 
the superconducting order parameter. Staggered super- 
conductivity does not induce magnetic order (either in- 
plane or perpendicular-to-plane). 

In summary, the Ginzburg-Landau theory of staggered 
superconductivity was applied to UPd2Al3 and UM2AI3. 
Motivated by odd-in-frequency pairing, we found the 
minima of the negative hopping model for Cooper pairs at 
the corners Hi,H 2 of the Brillouin zone. We introduced 
a two-component superconducting order parameter 771 , 772 
transforming as an ii-point irreducible representation of 
the space group. We constructed the free energy includ- 
ing gradient terms and coupling to symmetry-breaking 
strains. The superconducting transition is not split by 
any finite value of the strains. In the parameter space 
there are two zero- magnetic- field regions of stability char- 
acterized by the superconducting phases (771,7/2) = (i, 1) 
and (0,1), respectively. The upper critical field as a 
function of the temperature has no kink. We defined 
charge and magnetic order parameters and investigated 
whether the superconducting order parameter can induce 
them. For the phase (771,772) = (0, 1) there is an induced 



charge-density wave characterized by new Bragg peaks 
(■|tt,i, ini +712,77.3) which could be seen by neutron or 
x-ray scattering. 
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FIG. 1. Brillouin zone for the simple hexagonal Bravais 
lattice and elements of stars for representations characteriz- 
ing different order parameters. The stars {Mi,M 2 ,M 3 } and 
{H\,H2} correspond to the superconducting order parame- 
ters for UPt3 and UPd2Al3, respectively. The stars {F} and 
{Ki, K2} describe the charge and the magnetic order param- 
eters for UPd2Al3. 

FIG. 2. Original and cell-tripled hexagonal Bravais lattices 
(with only one plane shown) used to define the charge and 
the magnetic order parameters. The original lattice consists 
of all the points drawn. The new lattice consists of the points 
on the dotted lines and describes the most general charge 
and magnetic structures. Points characterized by different 
symbols represent uranium atoms with different charges and 
magnetic moments. 
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